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The phase diagram of a coupled tetrahedral Heisenberg model is obtained. The quantum 
chain has a local gauge symmetry and its eigenspectrum is obtained by the composition of the 
eigenspectra of spin-i XXZ chains with arbitrary distribution of spin-| impurities. The phase 
diagram is quite rich with an infinite number of phases with ferromagnetic, antiferromagnetic 
or ferrimagnetic order. In some cases the ground state and the low lying eigenlevels of 
the model can be exactly calculated since they coincide with the eigenlevels of the exactly 
integrable XXZ chain. The thermodynamical properties of the model at low temperatures 
is also studied through finite-size analysis. 

I. INTRODUCTION 

Recently an extensive amount of experimental and theoretical investigations on a large variety 
of quasi-one-dimensional mixed spin compounds have been reported in the literature l|- 3]. One 
of the main reason for the interest on these physical systems is due to the richness of their mag- 
netic properties, exhibiting antiferromagnetism, ferromagnetism, ferrimagnetism, magnetization 
plateaux, etc. 

The simplest model on this class is the mixed spin-(S', 5') Heisenberg model where we have 
two types of antiferromagnetically exchange-coupled spins S and S' located at alternated sites of 



a chain. According to the Lieb-Mattis theorem j^] in the cases where S ^ S' the model exhibits 
ferrimagnetic behavior, having a. {S — 5')-|-degenerated ground state, where L is the chain length. 
Extension of these models to quantum chains with two or more different quantum spins in the 
elementary cells increases considerably the number of quite interesting models, since the topology 
besides allowing a mechanism of frustration of spins also determines the ordering of the ground 
state, and the nature of the low-lying excitations. Several of these mixed-spin compounds exhibit 
local gauge invariance rendering a great amount of simplification on the analysis of these quantum 
chains through finite-size studies. In ^9] conformal invariance was exploited to obtain the critical 
properties of the simplest mixed-spin compound where we have at the odd sites a spin ^ connected 



with two disconnected spin i at even sites. In 



IC 



by introducing interactions among the spins 



2 



at the even sites it was shown that phase transitions occur as the strength of the interactions at the 
even sites changes. Another variant of these alternating compound spin chains is the orthogonal- 



dimer spin chain investigated in who is also known as the frustrated dimer-plaquete chain. All 
these quantum chains due to its gauge symmetry have a phase transition to a phase whose ground 
state in the thermodynamic limit is fragmented, i. e., it is formed by the composition of ground 
states of finite quantum chains with open ends. 

Motivated by these results we introduce and study in this paper a mixed-compound quantum 
chain that, as the previous models, exhibits a local gauge symmetry, but distinctly does not ex- 
hibit the fragmentation properties. These models describe the dynamics of spin-^ Pauli matrices 
located at the vertices of tetrahedrals and interacting as in Fig. 1. As we shall see this quantum 
chain reveals a quite rich phase diagram with ferromagnetic, antiferromagnetic, ferrimagnetic mas- 
sive phases as well massless disordered phases. These phases are connected through continuous 
and discontinuous phase transitions, as we change the coupling constants. The number of phase 
transitions can be finite or infinite, depending on the region of the phase diagram. 

Another interesting point about our model concerns the question of exact integrability of quan- 
tum chains. We show that part of the eigenenergies of our model are just the eigenenergies of 
th. exa^Untegrable spin-J XXZ ch.n. „h„« e.genene.gi« „e exactly .„„wn though .he Bethe 
ansatz |15l |. Depending on the exchange couplings these eigenenergies become the low lying ones, 
including the ground state energy. Consequently our model provides a clear example of a non- 
integrable model that contains part of its eigenspectrum exactly integrable. An interesting feature 
of our tetrahedral spin chain is the appearance of phases with a non-zero entropy per site in the 
thermodynamical limit, similarly as happen with the residual entropy of ice [l^ . 




FIG. 1: One dimensional array of tetrahedral spin cluster, where a and r are spin-i matrices, dashed lines 
correspond to the exchange interaction J in the triangle configuration and the solid lines represent the 
anisotropic interaction among the spins a and r spins, with anisotropy A. 



The paper is organized as follows. In §2 we present our model and by exploring its gauge 
symmetry we show that its eigenspectrum is given by the composition of the eigenspectrum of 
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spin-^ XXZ chains with spin-| impurities. In §3 we study the spectral properties of the spin-| 
XXZ chain in the presence of spin-| impurities. In §4 using numerical calculations and the results 
of previous sections we obtain the phase diagram of our tetrahedral spin chain introduced in this 
paper. In §5 we study the thermodynamical properties of the tetrahedral spin chain by evaluating 
its entropy, specific heat and the magnetic susceptibility at zero magnetic field. Finally in §6 we 
present our general conclusions. 

II. THE MODEL 

The model we introduce in this paper describes the dynamics of the spin-^ spins located in a 
one dimensional array of tetrahedral spin clusters as shown in Fig. ^ The spin operators attached 
to the basis and to the top of these tetrahedral are denoted by CTj {i = 1,2, . . . , ^) and rf, t\ and 
r? {i = 1,2, . . . , respectively. The interactions are among nearest neighbors in this geometrical 
structure (see Fig. P), and the Hamiltonian in a tetrahedral spin chain (L even) is given by 

L/2 
i=l 

where 

T, = {lT,Tr,T^)=Tt+T\ + Tl, (2) 



(<t„T,)a = afT^^ + a^ry + Aaf7;.^ (3) 

i = 1,2, . . . and we have imposed periodic boundary conditions. The spin operators {t^,t\,tI) 
have isotropic interactions among themselves (coupling J) and XXZ-like interaction (anisotropy 
A) with the spin operators a. The particular case where A = — 1 coincides with a special point of 
an alternated tetrahedral spin chain introduced in [l^ . The competitions between these two type 
of interactions will produce a rich phase diagram. Clearly if J is negative the T operators prefer 
to align forming a spin-| operator and we should expect an effective mixed anisotropic spin-(^, |) 
Heisenberg chain. However when J > we necessarily have frustrations among the r interactions 
in the basis of the tetrahedron, but the amount of these frustrations will depend upon the relative 
strength of the anisotropic interaction with the operators a. For \J\ < 1 we should expect that the 
XXZ interactions being stronger will induce, as in the case J < 0, an effective mixed spin-(|, |) 
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Heisenberg chain. On the other hand for J ^ 1 the interactions among the r operators, defined on 
the basis of the tetrahedron, are more important and in order to minimize their frustration they 
cluster into spin s = ^ operators and the Hamiltonian should behave effectively as an non-mixed 
and standard anisotropic spin-^ Heisenberg chain or XXZ quantum chain. Since the isotropic 
mixed spin-(|, ^) Heisenberg chain exhibits ferrimagnetic behavior in contrast with the standard 
antiferromagnetic behavior of the spin-^ XXZ chain, we should expect phase transitions lines 
separating the above regimes. 

The Hamiltonian exhibit a local gauge symmetry that will render a great simplification in 
our physical and numerical analysis of the model. It is simple to verify that the local operators 



G.=T2 = (T?+r^ + rf)^ i = l,2,...,-, (4) 

besides commuting among themselves commute with the Hamiltonian Consequently in the 
basis where all the {Gj} are diagonal the Hilbert space associated to can be separated into 
block disjoint sectors labeled by the eigenvalues {gi} of the operators {Gi} {i = 1,2, . . . , Since 
Ti is formed by the addition of three spin-^ angular momentum operators its eigenvalues are 
gi = S{S + 1), with 5 = |, ^. The eigenvectors corresponding to gi = in the T^-basis can be 
written as the quadruplet 



|3 3, 

,3 1, J_/ 



+ ++) 

+ + -) + I + + 



2'2/ 

I +) + + + 1 + 



,3 -l^ 1 / 



'2' 2 / V3V 

I- -) = I—) 
I2' 2 / I 

while the eigenvectors corresponding to = | can be written as the quadruplet 

I2' 2 / V2V 



(5a) 



fl +— ) - 1 — +)) (5b) 



and 



l5-l^)' = 7i(l + --> -2|- + -> + !-+». (50 
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As a consequence of the local gauge invariance the action of the Hamiltonian on the basis vectors 
only connects vectors with the same eigenvalue {gi}. Moreover in the case where 5i = | the first 
two eigenvectors given in H5b|) are not coupled through the interactions with the last ones, given 
in ((^ . This means that in fact the gauge choice = ^, for a given site i, can be chosen on two 
independent ways. 

Therefore using the basis © the Hilbert space associated to (0) are separated into sectors 
labeled by {51,92, • • • ,5^}, where gi = ^, |, |. For a given choice {gi} (z = 1, ...,-§) we verify 
that in this restricted subspace the effective Hamiltonian reduces to a general mixed anisotropic 
Heisenberg chain 



L/2 

} ~ ~ ^ {^2i-1^2i + ^2i^2i+l + '^2i-1^2i + ^2i^2i+l + ^{^21-1^21 + ^2i^2i+l) 
^ i=l 

i=l 

where at the odd sites we have spin-^ operators and at the even sites the spin operators Sj 

are spin-| or spin-^ Pauli matrices depending if = ^ or gi = |, respectively. Consequently a 

complete understanding of our model will demand the understanding of the infinite set of general 

mixed spin quantum chain ©. Particularly in order to calculate the phase diagram we should 

know, as a function of A and J, among all these 3^ mixed Heisenberg chains, which one provides 

the global ground state of the model. 

Some special gauge choices, according to our previous qualitative analysis are natural candidates 

for hosting the ground state in regions of the phase diagram (J, A). The gauge gi = g2 = • • • = 

gL = ^ give us the anisotropic mixed spin-(i, |) Heisenberg chain that at A = —1, due to the 
!U r 

Lieb Mattis theorem ]&\ exhibits a -g-fold degenerated ferrimagnetic ground state. Certainly for 
large and positive values of J this sector contains the ground state. 

The 2 2 distinct choices where gi = g2 = • • • = g L = i_i correspond to degenerated eigensectors 

2 ^ 

with the same eigenspectra as the standard and non-mixed spin-^ XXZ chain, whose exact solution 
is known '15] through the Bethe ansatz. 

It is interesting to stress here that although our coupled tetrahedral spin chain is not exact 
integrable it contains in its eigenspectrum an infinite set of exactly known eigenvalues. Moreover 
we should expect that this set of eigenenergies in special regions of the phase diagram (J, A) will 
contain the ground state energy and low lying excitations of the model (^. In these regions the 
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ground state of the model will be a 2 2 -degenerated ground state exactly known through the Bethe 
ansatz. 



Starting from the gauge gi 



(i = l,2 



; ^5 • • • ) 2 /' 



we may produce the other gauges by choosing 



some of gi = This gives an effective anisotropic spin-i Heisenberg chain with magnetic spin-| 



15 
4 • 



impurities, at the sites positions i where g., 

Differently from other quantum chains possessing local gauge symmetryj^ 



Id 
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1^ 



our Hamil- 



tonian 
models 
impurities. 



does not produce a fragmentation of the lattice. The fragmentation happens in those 
, because the effect of gauge fixing is to produce non-magnetic (spin zero) 




As we can see from @ the coupling J plays the role of a chemical potential since it gives the 
additional energy required to insert a spin-| impurity into the system. Before calculating the phase 
diagram of the Hamiltonian (jH)) we need to calculate, as a function of the impurity number, the 
ground state energy of the anisotropic spin-^ Heisenberg chain with spin-| impurities. 



III. THE HEISENBERG CHAIN WITH IMPURITIES 



The Heisenberg chain with impurities were studied in the case of periodic array of impurities 
by Fukui and Kawakami 18]. Their results reveal that whenever the two types of spins defining 
the chain are odd-half integers the system is gapless. In the present case, in order to calculate the 
phase diagram of (^, we need to know the lowest eigenenergy among all the ground-state energies 
of the spin-^ Hamiltonian with a given number of spin-| impurities located at the even sites. In 
general these quantum chains on a periodic lattice are not translational invariant. For example 
for L = 20 the configurations of impurities, apart from trivial rotations by two sites, are given 
by the configurations Cq, Ci, . . . , C4 in Fig. [2 Sind in the case of L = 16 the configurations with 
n = 4 impurities are given by Cq, Ci, ... ,Cj in Fig. |3| We calculate, using the Lanczos method the 
lowest eigenenergy of these Hamiltonians with arbitrary distribution of impurities, for lattice sizes 
up to L = 20 or up to L = 16 depending if the anisotropy A<— lor— 1<A<1, respectively. 
In the case of low density of impurities we are able to extend our calculations up to L = 24. 




FIG. 2: The configurations of two impurities of the chain with L = 20 sites. 






FIG. 3: The configurations of 4 impurities of the chain with L = 16 sites. 



The ground-state energy of the Hamiltonians with arbitrary distributions of n impurities (0 ^ 
n ^ ■^) is non degenerated for |A| < 1, two- fold degenerated for |A| > 1 and at A = —1, due to 
the Lieb-Mattis theorem |j8j|, the ground state is (2n + l)-fold degenerated. 

In order to compare the different ground-state energies for the Hamiltonians with distinct im- 
purity configurations we plot in Fig. 0J as a function of A, the energy differences Ec^ — Eqq 
{i = 1, . . . , 7) between the ground-state energies of the Hamiltonians with the impurity configura- 
tions of Fig. 01 for the chain will L = 16 sites. We see from this figure that for all values of A 
the ground-state energy Ec^, corresponding to the homogeneous distribution of impurities, has the 
lowest value, followed by the energies Ecf^ and Ec2- The energy Ecq, that corresponds to the con- 
figuration where the impurities are at the closest positions is the higher one. This behavior repeats 
for arbitrary lattice sizes and number of impurities, showing that due to quantum fluctuations the 
impurities repeal each other favoring the configuration where they are homogeneously distributed. 
Our numerical results also show that the minimum energy among these configurations E^{n, A, 0) 
decreases with the total number of impurities, which can be understood since the interaction be- 
tween two spins ^ and | give lower contributions to the energy than the interactions between two 
spins ^. 

The lowest eigenenergy E^in, A, J) of our tetrahedral Hamiltonian in the gauge sector with 



n values Qi = from ©, is then given by 

Elin,A,J) _ El{n, A,0) _ 3 1 _ 2n 
L ~ L 44 

From our previous discussions, these energies have a degeneracy 2 2"^" for |A| < 1, 2 • 2 2" 



(7) 



for 



A < -1, 2(^/2)2*-" for A > 1, (L + 1 + 2n){^l^) for A = 1 and (2n + 1) • 2*"" at A = -1. The 
ground-state energy of our tetrahedral spin chain for a given value of J and A, will be given 
by the lowest value of • In Fig. 13 we show as a function of J and at A = — 1 all the lowest 
eigenenergies El{n, —1, J)/L n = 0, 1, 2, . . . , for the lattice sizes L = 12, 18 and 20. They are 
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FIG. 4: The lowest energy of the Hamiltonian with L = 16 sites and impurities at the configurations C1-C7 
of Fig. 13 The energies are measured with respect to the lowest energy Eco of the Hamiltonian with the 
impurity configuration Cq. 



given by straight lines whose slope values depend upon the number n, labelling the gauge, and 
their intercept points with the vertical axis {E^{n, —1,0)/L) increase as the number of impurities 
decreases. Our numerical analysis shows that these plots depend strongly on the value of A. 

We see from Fig. |21 in agreement with our earlier physical arguments for A = —1, that 
for Ji(L,-l) < J < Ji(L, -1), with Ji(L, -1) -00 and Ji(L, -1) w 0.667, the ground 

2 4 2 4 

state energy is E^{L/2, —I, J) and belongs to the gauge sector gi = 92 = ■ ■ ■ = dL = 
that corresponds to the ground state of the anisotropic spin-(^, |) mixed Heisenberg chain. For 
Ji(L,— 1) < J < Ji(L,— 1), with Ji(L,— 1) 0.75, the ground state will be given by the gauge 

4 6 6 

sector with n = j values gi = and since among these gauge choices the lowest energy is given 
by the most homogeneous distribution of n = L/6 values ol gi = the ground-state energy 
is the same as that of the mixed spin-(^, ^, ^, |) Heisenberg Hamiltonian where we have spin-^ 
operators in all lattice points except at the equally spaced points i = 4j [j = 1,2,...). For 
Ji(L,— 1) < J < Ji(L,— 1) the ground state energy should come from the gauge sector with 

6 8 

the most homogeneous distribution with „ = f values of = i^. or equivalently it is given bv 

the mixed spin-(i, i, i, i, i, |) Heisenberg model. Our numerical results induce us to expect |l3i] 
that in fact for sufficiently large L, these crossings should continue, and we have in general for 
Jj_(L,— 1) < J < J 1 fL,— 1) (/ = 1,2,3,...) the ground-state energy in the sector with 
n = Yi values of i^j = or equivalently it is given by the ground-state energy of the mixed 
spin-(i^^' '^\^) Heisenberg chain. Fig. |Slalso shows that (J 1 (L,— 1) — Jj_(L,— 1)) ^ as / 

^ ^ 2(i + l) 21 

grows and for J > Jq{L, —1) the ground state energy is given by the gauge 51 = (72 = • • • = 5i = 
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that is equivalent to the ground-state energy of the exact solvable [l^ standard antiferromagnetic 
spin-^ Heisenberg chain. For other values of the anisotropy, as long A < 1, we have the same 
qualitative behavior as A = —1, with crossing values Jj_(L,A), however the crossing regions 

21 

5J^{L,A) = Jo(L,A) — Ji/2iL,A) decreases as we departure from A = —1. Particularly for 
A — oo or A > 1 all the crossings happen at the same point Jc(L, A) ( 6J'^{L,A) = 0), and in 
those cases the ground-state energy is the same as that of the anisotropic mixed spin-(^, |) Heisen- 
berg chain or the standard anisotropic spin-^ Heisenberg chain if J < Jc(L, A) or J > Jc(L, A), 
respectively. 



IV. THE PHASE DIAGRAM OF THE TETRAHEDRAL SPIN CHAIN 

In order to obtain the phase diagram of our Hamiltonian we need to estimate the bulk 
limit values Jr){oo, A) = limL^ooJr]{L, A) [rj = ^, j, ^, . . . ), where the ground-state change sectors 
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and phase transitions are expected. From the discussions of the last section and relation ((71) the 
crossing points for a finite lattice are given by 

^> - 3 ' 

with 

V' = (9) 

where r/' = |, g, |, . . . , 0. Consequently the phase transition lines Jr){oo, A) can be calculated from 
the finite size corrections of the energy levels i?^(r/L,A,0) of the mixed spin-(^,|) Heisenberg 
model with rjL spins | and (1 — r])L spins | homogeneously distributed. 

Based on general grounds the ground-state energy of these mixed Heisenberg models with a 
finite number n = rjL of spins | should behave asymptotically {L oo) as 

Elin,A,0)=e^Hm + FLin,A), (10) 

where 

Fi(n,A) = Foo(n,A) + o(l/L), (11) 

and e^^^(A) is the ground-state energy per site of the anisotropic spin-i XXZ chain, with 
anisotropy A, and Foo{n, A) plays the role of a surface energy (negative) due to the additional spin- 
I particles. In the region where the equivalent mixed Heisenberg chain is critical (— 1 ^ A ^ I'' 



due to the conformal invariance of the infinite system, the neglected term in pij) is 0{1/L) 
while in the non-critical regions it should behave as e~"^, where a > is proportional to the mass 
gap of the model. Moreover since the density of spin | is bounded, i. e., ^ r/ = ^ ^ i, we should 
expect in the limit where L ^ oo, n — > cxd but with V — x ^^P^ fixed that the surface energy per 
site is given by 

M^ = M,,A) „„,££(!^. (12) 

L/ L— >oo L/ 

where fooii], A) depends only on the value of rj and A. 

The finite-size extrapolations /oo(??, A) for the ferrimagnetic case (A = —1) for some densities 
are shown (diamonds) in Fig. ^ Although the curve formed by the union of these points is almost 
a straight line there exists a small decreasing upward concavity as r] decreases. This fact will be 
crucial to explain the existence of an infinite number of phase transitions for the model around 
A ~ — 1. In order to show these concavities let us consider the densities in the region j < ?7 < 5. 
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The largest lattice we can calculate with the most homogeneous distribution, for these densities is 
L = 16 (triangles) and L = 20 (circles). In Fig. Clby considering these lattice sizes we show the 
finite-size version of these derivatives, namely, the mass gap obtained by the increment of a single 
impurity in the system, i. e., E^{r]L, —1, 0) — E^{i]L — 1, —1, 0). This figure indicates that except 
for the densities close to rj = ^ ot rj = ^ this gap has the constant value Gi = —1.0019. The same 
behavior should also happens for the densities in other intervals 2(I+T) ^ ^ Ti ~ 2,3,...), 
where the gap has an almost constant value Gi that decreases as I increases. The small deviations 
shown in Fig. [7|for Gi around r] = 0.25 and ij = 0.5 are due to the finite-size effects. 




FIG. 6: Surface energy per site in the bulk limit fociv- ^1) a function of the density j] of spins | (top) 
and its derivative /^(fy, — 1) with respect to the density (bottom). 
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FIG. 7: Difference among the two lowest energies E^{riL, —1,0) and E'^ijjL — 1,-1,0) for the lattice sizes 
i = 16 (triangles) and i = 20 (squares) in the gauge sector where the density 77 = n/L of spins | changes 
from 0.20 to 0.5. The dashed line represents the bulk limit in the intermediate region -j < 77 < i. 



In the lower part of Fig. Elwe plot the numerical values obtained for Gi (Z = 1, 2, . . .). A similar 
analysis for other values of A show us that except for A > 1 (ferromagnetic region) a similar 
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behavior happens. These results imply that in the bulk limit we should expect that the surface 
energy per site is given by 

/oo(r/, A) = vGi{A), I = Int[l/2??]. (13) 

In the limit of lower density (rj ~ 0) and higher density (t? ~ ^) we can obtain, as a function 
of A, accurate results for Goo(A) and Gi(A), respectively. The results for ry ~ (r; ~ ^) are 
calculated from lattice sizes up to L = 24 {L = 20), where we have (L — 1) spins ^ (spins |) with 
a single impurity of spin | (spin ^). In Fig. IHlthe finite-size sequences are shown for the surface 
energy in the lower density limit rj = 0, and for some values of A. Using a cubic least-square fitting 
we obtained the extrapolated values of Goo (A) shown in tabled where the errors are estimated 
from the stability of the fitting when the largest lattices is neglected. A similar extrapolation for 
r/ ~ ^ shows that Gi(A) ~ —1.00187 for all values of A in the critical region. In the massive regions 
|A| > 1 we have two distinct behaviors, for A ^ — 1, Gi(A) = -|-A while for A > 1 (ferromagnetic 
region) Gi{A) = —A, for all values of I. This means that in the massive ferromagnetic regime 
/L(r/,A) = -r/A. 

The surface energy p3() when inserted in ((7| give us the bulk value (n — > oo,L — > oo, rj = n/L 
finite) for the energy per site of the lowest eigenenergy in the gauge sector with n = rjL values of 

e^iv, A, J) = lim^^oo ^"^^^^^"^^ = e^^^(A) + Q(A)r? - ^(1 - 2ry) J. (14) 

This last relation implies that the critical couplings J,,(cxd, A), given by ©, where the ground 
state changes its location on the gauge sectors are just given by 

J,(oo,A) = -^G_l(A), ^=1,1,1,... (15) 

In Fig. 1^ we plot for the densities rj = 0, |, ^ these lowest energies as a function of J for 
A = —1, A = and A = 1. We see from this figure that while for A = —1 and A = —0.5 there 
exist crossings from the curves with rj = and rj = ^, these do not happen for A > 1, in agreement 
with our discussions with finite size lattice (see Fig. It is important to stress here that those 
crossings are governed by the intercept with the vertical axis 600(7^, A), whose A-dependence in 
the bulk limit are given by the ground-state energy of the XXZ chain and the surface energies 
/oo(r/,A). The exact results of e^^^(A) |1^] and our numerical estimates for G;(A) show us the 
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FIG. 8: Extrapolation of the difference among the ground-state energy of the homogeneous anisotropic 
spin-^ Heisenberg chain and the spin-^ Heisenberg chain with a single impurity of spin |. 
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TABLE I: The energy gaps are tabulated at very low density (1 <C Z < oo) of concentration in the critical 
and massive antiferromagnetic regions. 

following asymptotic behavior. For A > 1 there is a single crossing from the density ry = ^ to = 
at A = where the ground state changes from that of the mixed spin-(^, |) Heisenberg chain 
to the standard spin-i XXZ chain. For A < 1 there exist level crossings for all the gauge choices 
with ?7 = 2j, ^ = 1,2,.... The range of values of J where these crossings happen is relatively small 
except around A « — 1. 
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FIG. 9: Phase diagram at zero temperature of the tetrahedral spin chain with Hamiltonian (Q. 

Collecting our numerical results we obtain the quite rich phase diagram shown in Fig. |^ for 
our tetrahedral spin chain. The phase diagram shows eight distinct phases that we are going to 
describe. The ground state energy in the phases Ai i, Di i,Fii (shadow region in Fig. |^ 

2 ' 2 2 ' 2 2 ' 2 

is given, apart from degeneracies, by the ground-state energy of the standard spin-^ XXZ chain, 
while the phases Ai a , Di s, Fi 3 are given by the ground-state of the anisotropic mixed spin-(i |) 

2'22'22'2 

Heisenberg chain. The intermediate phases Ai „ and Di „, actually are regions with an infinite 

2 ' ' 2 ' ' 

sequence of phase transitions where the ground-state energy is that of the anisotropic version of 
a mixed Heisenberg chain with spins ^ and |, where the density of spin | is rj = ^, ^, ^, . . . ,0. 
Despite of the large number of closed excited states on each of these phases the ground state is 
given by a distinct wave function with peculiar magnetic properties. The phase diagram of Fig. IHl 
can be divided in three distinct regions with respect to the variations of A. 

i) A ^ 1. For J > ^ (phase Fi 1) the system behaves as a ferromagnetic Heisenberg spin-| 
chain with z-magnetization = ^ ■§ while for J < ^ (Phase -Fi 3 ) it behaves effectively as 
a mixed spin-(^, |) Heisenberg chain in its two fully ordered state with total z-magnetization 
= ibL. At J = |A there is a first order phase transition separating these two phases. 

At A = 1 the phases Fi 1 and -Fi 3 end up into a gapless phase transition line where all the 

2 ' 2 2 ' 2 

lowest energies in the sectors with z-magnetization —M^ < < degenerate, rendering 
a ferromagnetic ordered state with total magnetization and a massless spectra. When 
A — > 00 (Ising limit) the ground-state energy per site in the phase Fi 1 is Eq = — 



2 ' 2 
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and in the phase -Fi 3 is Sq = —K- + 



3A I 3J 

4 8 



ii) |A| < 1. For J < | the system behaves as the disordered phase of the mixed spin-(^, |) 

Heisenberg chain. For J > Jo(oo,A) (shadow region) the model behaves as the standard 

and solvable XXZ chain in its massless regime. For i = Ji < J < Jo(oo, A) there exists an 

infinite number of phase transition separating the regions J^'(oo, A) < J < J,,(cx), A) (77 = ^, 

v' = 2(;+i) ' ^ ~ 2,3, ... ) where the model behaves effectively as a mixed spin-((i)^'~^, |) 

Heisenberg chain with density 77 = ^ (/ = 2, 3, . . . ) of spins |, in its massless regime. These 

intermediate phases happen in a small range of values of J in the phase diagram, except 

around A = — 1. At A = —1 all these phases end up into a phase transition line where the 

system try to establish an antiferromagnetic order. This antiferromagnetic order is achieved 

in the phase Di 1 . However in all the phases Di and Di 3 , due to the Lieb-Mattis theorem 
~ 2 ' 2 2 ' ' 2 ' 2 



the lowest eigenenergy in several sectors with distinct magnetizations (not all of them) 
are degenerated and we have phases with increasing ferrimagnetic order. 

in) A < — 1. For J < ^ (phase Ai 3) the model behaves as an effective spin (^, |) Heisenberg 

model with massive ferrimagnetic behavior, while for J > Jo{oo, A) (phase Ai 1 ) the model 

2 ' 2 

behaves as the antiferromagnetic spin-| Heisenberg model. For the intermediate values 
^ < J < Jriioo, A) the model exhibits a sequence of ferrimagnetic ordered phase, behaving 
effectively as a mixed spin-(^, |) Heisenberg model with distinct densities of spins |. 

It is interesting to stress here that in all the shadow region of Fig. ^ the ground state have a 
degeneracy of order 2 2 in a L-sites chain, due to the possible gauge choices of our model (|T|). This 
means that similarly as happen with the ice at low temperatures the entropy per site in all 
these regions is the finite number So = -f- In 2, where Kg is the Boltzmann constant. 

V. THERMODYNAMICS 

In this section we discus the thermodynamic properties of the tetrahedral spin chain The 
gauge symmetry of (0) enable us to write its partition function (see (fTHjl) 



i- 3. J 

^L(r,A,J) = (2e^'=B-)2^(__j ^ ZLiii,...,in;T,A), (16) 

n=0 l<ii<i2<---<in<L 

where Z^iii, . . . , i„; T, A) is the partition function of the spin-^ anisotropic Heisenberg chain with 
n spin-| impurities located at ii, . . . and Kg, as before, is the Boltzmann constant. 
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As T — > we can approximate 

g°(zi,...,in;A) 

ZL(ii,...,in;T,A)c^gl{H,...,in)e , (17) 

where E^{ii, . . . , A) is the ground-state energy of the equivalent spin-^ Heisenberg chain with 
impurities at {ii, . . . , whose degeneracy is given by gi{hi ■ ■ ■ , in)- Even with the approximation 
H17|) the evaluation of Zl{T) is not simple for general values of J and A. In the ferromagnetic 
regime A > 1 however some additional simplifications occur. As discussed in §2 the degeneracy 
g^{ii, . . . ,in) = 2 and the ground-state energy . . . , in,; A) does not depend on the specific 

locations of the impurities, only on their number. Moreover as L — > oo, taking into account that 
for A > 1, G/(A) = - A in we have 

Elih, ...,in;A) = E\n, A) = Le^^^(A) - nA. (18) 

Inserting ((TH|) and ((T7|) in (fTB|) we obtain 

ZL(r,A,J)=2e "s"^ 2e*''fl^ +e(^+4-^)/''s^ % (19) 

for T — > 0, -L ^ CO and A > 1. This expression gives the entropy per particle 

3J-2A 

s{T, A, J) = ^ ln(2 + e-"^) + {3J_-2A) _ ^^^^ 



2 + e 



T 



r = 0.001/k^ 




FIG. 10: Entropy at low temperatures in the ferromagnetic regime (A = 1, L — > oo) and in the ferrimagnetic 
regime (A = -1,L = 12). 



Exactly at T = we have s{T,A,J) = for J < jj} = |A and s{T,A,J) = ^ln2 for 
J > Jo{oo, A). The expression (|^n|) should be valid for k^T < 0.1, since the gap in the ferromagnetic 
regime, that controls the approximation (jlZj) . is of order of 0.1 for the XXZ chain with an arbitrary 
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number of impurities. In figure ITUl we show the entropy per site given by H20() at T = 0.001/tg and 
A = 1. We see from this figure that the net effect of the temperature is to produce a peak at the 
zero temperature quantum critical point Jo{oo,A). 

The entropy for the other anisotropics A < 1 can only be estimated numerically from the 
exact numerical calculations on small lattices. We show in figure ll()( in order to compare with 
the ferromagnetic case A = 1, given by (|20j) . the numerical results for the lattice L = 12 of the 
ferrimagnetic model at A = — 1 and T = 0.001/ Kg. In this case the entropy shows small peaks 
around the zero temperature phase transition couplings Ji (oo, A), Ji (oo, A), .... Strictly at T = 

2 4 

and for —1 < A < 1 these peaks disappear and the entropy per site is given by the set of plateaux 



s(0, A, J) = ^(1 + -)ln2, Ji (oo.A) < J < J i foo. A). (21) 
As we can see in the phase diagram of Fig. |^for A < 1 the sizes of the plateaux (J i (oo, A) — 

2(i + l) 

J_i_ (oo, A)) — > 0, for any value of A < 1, while for A > 1 and A <C — 1 there exists a single plateaux 
starting at J = Jo(oo, A) = |A. 

The evaluation of the partition function (jlZj) for arbitrary temperatures can be done by col- 
lecting all the eigenvalues of the several XXZ quantum chains with arbitrary distribution of spin-| 
impurities at the even sites of the lattice. We can perform these calculations up to the lattice size 
L = 10. In order to extend these calculations for L = 12 we would need to diagonalize matrices of 
order 32424 x 32424, which is beyond the computational possibilities nowadays. 



A. The specific heat 

Let us consider initially the specific heat at very low temperature. In the massive ferromagnetic 
case we have from ()2U() 

/'\T—9A\2r 3J-2A 3J-2A 

Cl{T,A,J)=kJ ''' ) 2e"^«s-+e-B- 

for L oo, KgT < 0.1 and A > 1. In order to have an idea of the finite size-effects in the low 
temperature regime, where they are more important, we plot in Fig. ^2 the bulk limit (L oo) 
of the specific heat at A = 1 and k^T = 0.01, given by (|22)) (solid line) together with its finite 
size versions for lattice sizes L = 6, 8 and 10 (dotted, dashed and dotted-dashed lines respectively), 
obtained by a direct numerical diagonalization of the quantum chains. These curves show that 
the main features of the specific heat for the ferromagnetic chains at T = 0.01 are already shown 



(22) 
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for relatively small lattices. From l|22|) for L — > cxo, A > 1 and k^T < 0.1 two peaks happen at 
J = ji^) = |A-1.769772K3rand J = J(2) = |A + 1.485203Kgr with intensities Coo (T, A, = 
0.388090112«;^ and Coo(T, A, = 0.12038896k^ . It is interesting to observe that for these low 
temperatures the intensity of these peaks does not depend on the temperature and they happen 
nearly at the zero temperature phase transition point Jo(cxd, A) = |A. Also we can see from 
that the peaks at J^^^ and J*-^-* are due to the temperatures fluctuations of the ferromagnetic XXZ 
chain levels and the ferromagnetic mixed spin chain, respectively. Strictly at T = these two 
peaks add up at the phase transition point J = Jo(oo, A) = |A. 

In the case of the non- ferromagnetic regime (A < 1), even in the low temperature region, an 
analytical calculation of the partition function, or the specific heat, is quite complicated. In order 
to compare with the ferromagnetic case we plot in Fig. ^Jthe specific heat of the ferrimagnetic 
chain A = — 1 at k^T = 0.0005 as a function of A, for the lattice size L = 12. This calculation 
was performed directly from the partition function H16|) obtained by considering only the lowest 
eigenenergy of the several eigensectors. This lattice size give us a good finite-size realization for the 
densities V — \ ^-^d g, and comparing with similar results for L = 16 (not plotted in Fig. Ill|l we 
verify that the main features of the L — > oo system is already kept with this lattice size. We see 
that as in the ferromagnetic case the peaks happen around the zero temperature phase transition 
points (see the phase diagram of Fig. 




J J 

FIG. 11: The specific heat at very low temperatures for the finite rings L — 6,8,10 in the ferrimagnetic 
(A = —1) region and in the ferromagnetic region (A = 1) with L = 10 sites. 

In order to study the specific heat properties at higher temperatures we diagonalize numerically 
the full eigenspectrum of the quantum chain with L = 6, 8 and 10, and compute their specific 
heat. These results are shown in Fig. El for the ferrimagnetic models A = — 1, with the values of 
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J = |, I and — |, and for A = —10 for J = 5.5 and 7.5. These curves except at low temperatures 
give us a good idea about their convergence to the bulk limit L oo. At low temperatures these 
curves display some anomalous behavior, with small peaks and inflexion points. These anomalous 
behavior are due to the combination of the effect of the zero temperature phase transitions, that 
is not quite sensitive to the lattice size, with the normal finite-size effects, as normally happens 
when we study a finite chain that becomes massless as L — > cxd. These effects can be distinguished 
(see Fig. [T^ as the lattice size grows (see also Fig [TT|) . We also show in Fig. [T^some Cl(T, A, J) 
for the ferrimagnetic point A = —10. In this case as we can see in Fig. El we have approximately 
a direct first order phase transition happening at J = Jo(cxd,A) = — |A = Consequently 
for J = Jo(oo,A) we should observe the anomalous finite-size behavior as mentioned before (see 
the small peak at J ~ 5.5 and T ~ 0.02), that disappears as J increases. In the vicinity of 
J = Jo(cxD, A) we observe a small peak which is a signature of ferrimagnetism. As we departure 
from the phase transition region J ^ Jo(cxd, A) this small peak disappears and we have only the 
standard Schottky-like peak of the specific heat. 
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FIG. 12: The behavior of the specific heat for the finite rings L — 6,8, 10 in the ferrimagnetic (A 
region and in the massive antiferromagnetic region (A — —10). 



-1) 



From now on we restrict our calculations with the finite ring with L = 10 sites. As illustration 
we plot in Fig. ^Jthe specific heat as a function of J for some temperatures and for the anisotropies 
A = — 1 and A = —2. For other values of A the general feature of two maximums repeats. The 
reason for the occurrence of these two maximums can be understood from Fig. 1111 The sharp peak 
that happens in Fig. ^2 for low temperatures is smoothed as the temperature increases. Similarly 
as in Fig. ^Jthe minimum between the two peaks happens in the region where we have the series 
of phase transitions for the model at zero temperature, i. e., 0.6667 < J < 0.8063 for A = —1 and 
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1.333 < J < 1.369 for A = —2 (see Fig. EJ. The maximums appearing for J < and J > Jc 
are due to the eigenlevels of the spin-i XXZ and spin-(i, |) Heisenberg chain, respectively. Both 
maximums diminish as |J| ^ oo recovering the standard specific heat of the spin-| Heisenberg 
chain (J —>■ oo) or spin-(i, |) Heisenberg chain (J — oo). 

A = -l A = -2 




FIG. 13: Specific heat as a function of J for fixed temperatures in the ferrimagnetic (A = —1) and antifer- 
romagnetic (A = —2) regions. 



In order to give an overview about the specific heat in both regions of ferromagnetic and 
ferrimagnetic behavior we plot in Fig. the specific heat as a function of the temperature for 
several values of J and for the anisotropy values A = — 1 and A = 1. At the limits J ^ oo 
and J — oo the specific heat coincides with that of the spin-| Heisenberg chain and mixed 
spin-(i, |) Heisenberg chain, respectively. Let us consider the ferrimagnetic case A = — 1. For 
J = 10 ( J = —10) we see for k^T < 1.5 the specific heat behavior of the spin-i Heisenberg chain 
(spin-(|, |) Heisenberg chain), whereas for k^T > 1.5 it appears a peak with maximum k^T ~ 6. 
For J = 5 ( J = —5) and k^T < 0.5 the behavior is that of a spin-i Heisenberg chain (spin-(i, |) 
Heisenberg chain), whereas for k^T > 0.5 the curves deviate from the J — > oo (J ^ — oo) cases 
to form another peak with maximum at /tgT ~ 2.3. As we can see in Fig. E]a similar behavior 
appears for the ferromagnetic cases where A = 1. The second maximum appearing in Fig. 1141 can 
be understood from the limiting case of the model where \J\ ^ A. In this case we should consider 
only the interactions of the spins around the triangles of Fig. ^ {t^,t^, Tf) and we have a set of 
decoupled 3-sites chain whose specific heat is given by 

/ 8 7 \ 2 / 3J -3J X -2 

Coo(r,A,J)=2K,(^) (e-s-+e-B-) . (23) 

This relation gives a good approximation for the second maximum appearing for positive and 
negative values of J. If we compare this last relation with H22|l , which is valid for the ferromagnetic 
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model at low temperatures we verify that the fa 
the presence of two distinct maximums, now in 
the decoupled system an irrelevant contribution. 




FIG. 14: Specific heat as a function of T for fixed 
and in tlic ferromagnetic (A — 1) cases. 



ctor 2 appearing in H22() that was responsible for 
((2ni) is missing since the degeneracy 2 2" gives for 




2 4 ^6 8 10 



of the parameter J in the ferrimagnetic (A = — 1) 




In Fig. El we illustrate the distinct behavior of the specific heat, as a function of temperature, 
in the several phases of the model. The curves are for the L = 10 chain and are plotted for a fixed 
value of the parameter J = 0.75 and for some values of the anisotropy A. All curves show beyond 
the standard high temperature broad peak another one at low temperatures due to the proximity 
of the zero temperature phase transitions. The hight of these peaks depend on their location relate 
to the nearby zero temperature phase transition points (see FigEland FiglTTjl. 
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B. The magnetic susceptibility 

In order to calculate the magnetic susceptibility of our tetrahedral spin chain we introduce 
in the additional term 

L/2 

B = -5/x^/i^(af + (24) 

i=l 

where h is an external magnetic field, g is the Lande factor and jig is the Born magneton. This is 
equivalent to the addition of the term 

L/2 

B = -5/i,/ij;(cT2Vi + S|,) (25) 

i=l 

in the equivalent mixed spin Hamiltonian (jEJ. 

We calculate in this section the zero-field magnetic susceptibility per site xl{T., J) for our 
tetrahedral spin chain, i. e., 

Xl(T, a, J) = '^^InZUT, A, J, h)k=o, (26) 

where 2l{T, A, J,h) is the partition function with the inclusion of the magnetic interactions 
(PI or (Uni)- 

Let us consider initially the low and high temperature limits for the magnetic susceptibility. At 
low temperatures the degeneracy d{J, A) of the ground state, that depends on the value of J and 
A (see the phase diagram in Fig. |^ give us the behavior 



lunrx,(r,A,J) = ^!^^^|±^, d{J,A) = 2m + l. (27) 

On the other hand at high temperatures xl(T, A, J) should be independent of the field strength 
and saturate the Curie law value for independent spins. The magnitude of the independent spins 
depend on the value of the parameter J. In the cases where | J| ^ A we can neglect the anisotropy 
A in © and we have the asymptotic results 



TxLiT, J) = ^(1 + —^^ —). (28) 

e 



- *«B + e B ^ 



This means that when J — > — oo we recover the Curie law for the spin-(i, |) Heisenberg mixed 
chain: Txl{T, A, J) = while for J — > oo we obtain the Curie law for the spin-i XXZ 
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chain: Txl{T,A,J) = At finite values of |J|, but T ^ cx), we obtain from ^ a third 

^ 2 2 

possibiUty: Txl{T, A, J) = ^ ^ . We remark that the values of the constants in these Curie laws 
are independent on the existence of an infinite number of phase transitions or a single one, as it 
happens in the ferromagnetic regime. 

We calculated numerically the magnetic properties of our system through a direct diagonaliza- 
tion of the Hamiltonian for finite chains up to L = 10 sites. In Fig. we show Txl{T, A, J) 
as a function of the temperature for the finite chains with L = 6, 8 and 10 sites in the ferrimag- 
netic case (A = —1) and the values of J = 1, | and 0. At J = 1 the product Txl{T, A,J)—^0 
when T ^ 0. In this limit it behaves as the spin-^ Heisenberg chain whose ground-state energy 
is non-degenerated (m = 0), in agreement with the relation (|27() . On the other hand at J = 2/3 
we have a first order phase transition at T = connecting the phases with similar behavior as 
the spin-(i, ^, ^, |) and spin-(i, |) Heisenberg mixed chain. Finally at J = we are in a regime 
with similar behavior as the mixed spin-(^, |) Heisenberg chain with a ground-state degeneracy 
d{A, J) = d{—l,0) =2^ + 1. In this case, at low temperatures, we obtain the limiting value 

2 2 

TxUT, A, J) ^'^-^ predicted by (|17|), that in the bulk hmit goes to infinity. 

A=-l 

1 

0.8 
5 0.6 

^-0.4 
0.2 

^0 0.1 0.2 0.3 0.4 0.5 

T 

FIG. 16: The magnetic susceptibility as a function of the temperature for finite rings in the ferrimagnetic 
(A = —1) region. 

From now on we discuss the magnetic properties only for the finite-size chain with L = 10 sites. 
In Fig. El we plot Txl{T, A, J) as a function of the temperature for some values of J, in the 
ferrimagnetic case (A = —1) and in the ferromagnetic case (A = 1) with zero external magnetic 
field. Let us consider initially the ferrimagnetic case (A = —1). The limiting case J — > oo (long- 
dashed line) give us the spin-(^,|) Heisenberg chain with a ground-state degeneracy d{A,J) = 

1,— oo) = 2^ + 1. As the temperatures increases the magnetic susceptibility decreases with 
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a minimum value at n^T ~ 1.0 and then increases towards the value | — j^, in agreement with 
H28|). At J = (doted-dashed line) we have a similar behavior as the spin-(^, |) Heisenberg chain 
at very low temperatures whereas for high temperatures the system behaves, in agreement with 
(^5]) . as Txl{T, A, J) — > ^ ^ ^ (finite J). We also show in this figure the curve with J = | 
(solid line), i. e., the zero temperature transition point connecting the phases Ji and Ji. For 

2 4 

J = 0.8 (doted-doted-dashed line) we see in the figure that the low temperature behavior tends 
toward the behavior of the J — > — oo case, i. e., the standard behavior of the XXZ chain, where 
Txl{T, a, J) — > in agreement with (|27l) . since the degeneracy of the ground state is now finite. 

In order to compare with the ferromagnetic chains we also show these curves for the ferro- 
magnetic case. The degeneracy at T = is now large (infinite for L oo), but similarly as in 
the A = — 1 case depends on the value of A (see the phase diagram of Fig. The magnetic 
susceptibility in the ferromagnetic regime decreases monotonically at high temperatures tending 
towards the Curie law behavior with constants given by (|28|1 . 

A = -l A==l 




%' I 2 3 4 5 °0 0.5 L 1.5 2 

FIG. 17: The magnetic susceptibility as a function of the temperature for fixed values of J, in the ferrimag- 
netic (A —1) and ferromagnetic (A = 1) regions. 



VI. CONCLUSIONS 

The phase diagram of the coupled tetrahedral Heisenberg chain given by the Hamiltonian Q 
is obtained. Exploring the local gauge symmetry of the model we show that its eigenenergies is 
given by the addition of the eigenenergies of the XXZ chain with an arbitrary distribution of spin-| 
impurities at the even sites of the lattice. This imply that part of its eigenspectra is just given by 
the eigenenergies of the exactly integrable XXZ chain. Moreover, in some phases the ground state 
energy as well the low lying energies coincides with the eigenlevels of the standard XXZ chain. 
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and we have an example of a non-integrable model where the ground-state, as well the low-lying 
eigenenergies are exactly known (there are an infinite number of them in the thermodynamic limit). 

The pure XXZ chain exhibits a massless disordered phase for — 1 < A < 1 and massive phases 
ferromagnetically (A > 1) or antiferromagnetically (A < —1) ordered. The addition of spin-| 
impurities on the XXZ chain brings the possibility of ferrimagnetic order on the system due to the 
Lieb-Mattis theorem jS] • As a consequence the phase diagram of the Hamiltonian (see Fig. is 
quite rich, exhibiting an infinite number of phase transitions connecting phases with all the above 
types of long range magnetic order. In the intermediate phases the ferrimagnetic order increases 
until the saturated antiferromagnetic order is reached. It is interesting to note that in the phases 
^1 i) Di ir,, Fi 1 as well Di (see Fig. El the degeneracy of the ground state is extensive with 

2 ' 2 2 ' , 2 ' 2 2 ' ' 

the lattice size. These degenerated eigenenergies are connected even at zero temperature through 
quantum fluctuations and we should expect co nseq uently a finite residual entropy per site at zero 
temperature, similarly as happens with the ice 

The thermodynamical properties of our tetrahedral spin chain was also studied at low temper- 
atures. Thanks to the local gauge symmetry of the model the partition function at an arbitrary 
temperature was obtained exactly for chains with up to 10 sites. The specific heat as a function of 
the temperature in general exhibits besides an standard Schottky-like peak an additional peak, due 
to the additional eigenspectrum contribution to the density of states coming from the distinct gauge 
sectors. Close to the quantum phase transition couplings these peaks fuse into a single one. The 
magnetic susceptibility at zero magnetic field and low temperatures also presents some interesting 
features. As we change the exchange parameter J it shows a minimum for the antiferromagnetic 
case and a broad maximum for the ferromagnetic case. Moreover at high temperatures, depending 
on the value of J, we obtain the Curie law behavior with three distinct constants. 
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